Renewable  and  Sustainable  Energy  Reviews  33  (2014)  636-644 


ELSEVIER 


Contents  lists  available  at  ScienceDirect 

Renewable  and  Sustainable  Energy  Reviews 

journal  homepage:  www.elsevier.com/locate/rser 


Heat  transfer  enhancement  in  parabolic  trough  collector  tube  using 
Al203/synthetic  oil  nanofluid 

T.  Sokhansefat* 1 2 3 4 5 6,  A.B.  Kasaeian3’*  F.  Kowsary 

a  Department  of  New  Sciences  and  Technologies,  University  of  Tehran,  Tehran,  Iran 
b  School  of  Mechanical  Engineering,  University  of  Tehran,  Tehran,  Iran 


CrossMark 


ARTICLE  INFO 


ABSTRACT 


Article  history: 

Received  21  September  2013 
Received  in  revised  form 
14  January  2014 
Accepted  14  February  2014 
Available  online  12  March  2014 

Keywords: 

Parabolic  trough  collector 
Non-uniform  heat  flux 
Monte  Carlo 
Nanofluid 


Three-dimensional  fully  developed  turbulent  mixed  convection  heat  transfer  of  Al203/synthetic  oil 
nanofluid  in  a  trough  collector  tube  with  a  non-uniform  heat  flux  was  numerically  studied.  The  effect  of 
AL203  particle  concentration  in  the  synthetic  oil  on  the  rate  of  heat  transfer  from  the  absorber  tube  was 
also  investigated.  Various  nanoparticle  concentrations  (  <  5%  in  volume)  at  the  operational  temperatures 
of  300  K,  400  K,  and  500  K  were  used  in  the  current  study.  The  heat  flux  distribution  on  the  outer  surface 
of  the  absorber  tube  was  non-uniform  in  the  circumferential  direction  but  uniform  in  the  axial  direction. 
The  heat  flux  in  the  circumferential  direction  was  obtained  using  the  Monte  Carlo  ray  tracing  technique. 
Three-dimensional  Navier-Stokes  mass,  momentum  and  energy  equations  were  solved  using  the  FLUENT 
software.  A  second  order  discretization  was  used  for  the  convective  and  diffusive  terms  and  the  SIMPLE 
scheme  was  applied  for  the  velocity-pressure  coupling.  The  outlet  temperature  of  the  absorber  tube 
from  the  current  numerical  simulations  was  compared  with  existing  experimental  data.  The  absolute 
errors  were  found  to  be  less  than  3.8  C.  The  current  numerical  results  show  that  the  convection  heat 
transfer  coefficient  has  a  direct  dependency  on  the  volumetric  concentration  of  nanoparticles  in  the  fluid. 
In  addition,  the  heat  transfer  enhancement  due  to  the  nanoparticles  in  the  fluid  reduces  as  the  absorber 
operational  temperature  is  increased. 
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1.  Introduction 

The  global  energy  demand  is  increasing  in  an  accelerating  rate. 
Fossil  fuels  (oil,  coal  and  natural  gas)  are  currently  supplying  over 
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90%  of  the  world's  energy  demands.  The  fossil  fuels  are  non¬ 
renewable  and  limited  in  supplying  energy  indefinitely.  These 
energy  resources  will  eventually  be  depleted.  In  addition,  combus¬ 
tion  of  fossil  fuels  makes  carbon  dioxide  emission,  which  is  one  of 
the  greenhouse  gases  contributing  to  global  warming.  Use  of 
renewable  energy  resources,  on  the  other  hand,  can  lead  to  reduction 
in  fossil  fuel  consumption  and  in  turn  protection  of  the  environment. 
Solar  energy  is  one  of  the  largest  renewable  energy  resources. 
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Nomenclature 

A  T 

Average  fluid  temperature  above  ambient  air  tempera- 

ture(K) 

cP 

Specific  heat  capacity  (J/kg  K) 

u,v,w 

Dimensionless  velocity  (m/s) 

g 

Gravitational  acceleration 

9 

Angular  coordinate 

h 

Heat  transfer  coefficient  (W/m2  K) 

r 

Radial  direction 

I< 

incident  angle  modifier 

A 

Dynamic  viscosity  (kg/m  s) 

K 

Thermal  conductivity  (W/m  K) 

4> 

Particle  volume  fraction 

L 

Length  of  tube  (m) 

p 

Density  (kg/m3) 

P 

Pressure  (Pa) 

a 

Absorber 

Q" 

Heat  flux  (W/m2) 

bf 

Base  fluid 

Da 

Diameter  of  absorber  tube  (m) 

/ 

Focal  length  (m) 

Re 

Reynolds  number 

eff 

Effective 

T 

Temperature  (I<) 

s 

Solid  particle 

w 

Wall 

The  greatest  advantage  of  the  solar  energy  over  alternative  forms  of 
fossil  energy  resources  is  that  it  is  clean  and  abundant. 

Although  the  costs  of  solar  energy  have  gone  down  and 
continue  to  fall,  the  levelized  costs  of  solar  energy  are  still  much 
higher  than  conventional  energy.  The  levelized  cost  of  solar 
concentrated  power  (CSP)  ( ~  194-336  USD/MW  h)  is  four  times 
than  that  of  supercritical  coal  without  carbon  capture  and  storage 
(~43-88  USD/MW  h)  [1],  Cost  is  one  of  the  major  factors  inhibit¬ 
ing  development  of  trough  collectors.  Increasing  the  efficiency  of 
the  solar  collector  systems  can  in  part  alleviate  this  problem.  Many 
attempts  have  been  made  to  improve  the  performance  of  these 
systems  [2-5]. 

Heat  transfer  can  be  enhanced  by  increasing  the  thermal 
conductivity  of  the  heat  transfer  fluid  (HTF).  Thermal  conductivity 
of  metallic  particles,  metallic  oxides  and  nanotubes  is  relatively 
higher  than  that  of  liquids.  Addition  of  fine  particles  into  heat 
transfer  fluids  (thus  forming  nanofluids)  can  significantly  increase 
the  heat  transfer  rate  [6-8],  It  is  obvious  that  the  researches  on 
the  applications  of  nanofluids  have  been  popularized  during  the 
recent  years.  Many  studies  have  been  carried  out  on  the  effect  of 
nanofluids  and  nanofluid  size  on  heat  transfer  [9-10],  A  number  of 
investigations  show  the  improvement  of  convection  heat  transfer 
using  nanofluid.  Choi  et  al.  [11]  studied  engine  oil/carbon  nano¬ 
tube  nanofluids.  They  observed  160%  enhancement  of  thermal 
conductivity  using  1%  volume  carbon  nanotube.  The  concentration 
of  nanoparticles  has  high  effects  on  the  fluid  convection  heat 
transfer  coefficient  [12],  Akbari  and  Behzadmehr  [13]  investigated 
the  mixed  convection  heat  transfer  of  a  water/Al203  nanofluid 
system  in  a  horizontal  tube  with  uniform  heat  flux.  Their  results 
demonstrate  that  the  convective  heat  transfer  coefficient  is 
increased  from  9%  to  15%  using  2%  and  4%  volume  A1203  nano¬ 
particles,  respectively.  Duangthongsuk  and  Wongwises  [14]  stu¬ 
died  a  system  of  water  nanofluid  consisting  0.2%  volume  Ti02 
nanoparticles  which  caused  6-11%  enhancement  of  heat  transfer 
coefficient. 

Teng  et  al.  [15]  surveyed  the  changes  in  heat  transfer  of  A1203- 
water  nanofluid  at  different  diameter  sizes  of  nanoparticles.  They 
declared  better  thermal  conductivity  in  the  smaller  nanoparticle 
diameters.  In  another  study,  Heris  et  al.  [16]  compared  the  heat 
transfer  performance  of  Al203/water  with  CuO/water  nanofluids. 
They  observed  that  the  obtained  heat  transfer  enhancement  with 
A1203  nanoparticles  is  higher  than  that  of  CuO  nanoparticles.  As  an 
example,  they  stated  that  at  a  Peclet  number  of  5000,  the  heat 
transfer  coefficient  enhancement  ratio  was  1.29  for  2.5%  volume  of 
A1203  nanoparticles  whereas  it  was  1.23  for  CuO/water  nanofluid 
for  the  same  particle  volume  fraction.  Bianco  et  al.  [17]  considered 
a  laminar  flow  of  Al203/water  nanofluid.  Under  constant  wall  heat 
flux  boundary  conditions,  they  analyzed  the  heat  transfer  of 


nanofluid  using  both  single  phase  and  two-phase  approaches. 
They  noted  that  the  difference  between  the  results  of  single  phase 
and  two-phase  approaches  is  small,  especially  when  the  tempera¬ 
ture  dependency  of  the  thermo-physical  properties  is  taken  into 
account. 

Using  nanofluid  as  a  working  fluid  in  solar  systems  is  a  novel 
approach  to  increase  the  efficiency  of  solar  systems.  Lenert  and 
Wang  [18]  studied  the  influence  of  different  variations  in  nano¬ 
fluid  volumetric  receivers  both  theoretically  and  experimentally. 
In  their  theoretical  part,  a  one-dimensional  transient  heat  transfer 
model  was  supposed  and  the  enhancement  of  receiver  efficiency 
with  augmentation  of  nanofluid  height  (H)  and  incident  solar 
flux  was  proved.  In  the  experimental  part,  carbon  coated  cobalt 
nanoparticles  were  added  to  Therminol®  VP-1  in  a  liquid-based 
volumetric  receiver.  They  illustrated  that  the  efficiency  of  nano¬ 
fluid  volumetric  receivers  is  increased  with  increasing  solar 
concentration  and  nanofluid  height.  Also,  they  stated  that  the 
receiver-side  efficiencies  are  predicted  to  exceed  35%  when  nano¬ 
fluid  volumetric  receivers  are  coupled  to  a  power  cycle  and 
optimized  with  respect  to  the  optical  thickness  and  solar  exposure 
time.  Lu  et  al.  [19]  evaluated  the  influence  of  water-based  CuO 
nanofluids  upon  an  open  thermo  syphon  utilizing  in  a  high- 
temperature  evacuated  tabular  solar  collector  (HTC).  Compared 
to  water,  nanofluid  improved  the  thermal  performance  of  the 
evaporator;  also  a  30%  enhancement  was  observed  for  the  evapor¬ 
ating  heat  transfer  coefficient. 

Yousefi  et  al.  [20]  investigated  the  effect  of  A1203-H20  nano¬ 
fluid  on  the  efficiency  of  a  flat-plate  solar  collector  experimentally. 
In  comparison  with  pure  water,  they  declared  28.3%  enhancement 
of  efficiency  for  0.2  wt%  nanoparticle  concentration.  Also  they 
showed  that  utilizing  MWCNT-H20  nanofluid  enhanced  the  effi¬ 
ciency  of  a  flat-plate  solar  collector  [21],  Khullar  et  al.  [22]  studied 
the  nanofluid  usage  in  a  concentrating  parabolic  solar  collector 
theoretically  and  compared  the  results  with  the  experimental  data 
of  the  conventional  concentrating  parabolic  solar  collectors.  Their 
results  demonstrate  5-10%  higher  efficiency  as  compared  to  the 
conventional  models.  Kasaeian  et  al.  [23]  investigated  the  heat 
transfer  enhancement  for  Al203/synthetic  oil  nanofluid  in  a  para¬ 
bolic  trough  collector  tube  numerically.  Their  results  show  that  the 
heat  transfer  coefficient  is  increased  as  the  concentration  of  the 
nanoparticles  in  the  base  fluid  is  increased.  Otanicar  et  al.  [24] 
examined  the  effect  of  different  nanofluids  (carbon  nanotubes, 
graphite  and  silver)  on  the  performance  of  a  direct  absorption 
collector.  The  efficiency  rose  with  the  enhancement  of  particles 
concentration  but,  after  a  volume  fraction  of  5%,  the  efficiency 
diminished  slightly.  Since  the  nanofluid  receiver  can  potentially  be 
either  useful  in  a  direct  steam  generation  system  (DSC),  some 
works  are  reported  in  the  literature  [25]. 
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Table  1 

Economic  comparisons  for  conventional  and  nanofluid-based  solar  collectors  [26], 


Parameter  Conventional  solar  Nanofluid  solar 

collector  (S)  collector  (S) 


Capital  cost 

Independent  costs 

200 

200 

Area  based  costs 

397.8 

327.8 

Nanoparticles 

188.79 

Total  capital 

597.8 

716.59 

(one  time  cost) 

Total  maintenance 

96.23 

115.35 

(for  15  year  life) 

Total  cost 

694.03 

831.94 

Electricity  cost 

270.13 

278.95 

saving  per  year 

Years  until  electricity 

2.57 

2.98 

saving = cost 

Natural  gas  cost 

80.37 

83.02 

saving  per  year 

Years  until  natural 

8.64 

10.02 

gas  saving = cost 

Electricity  price 
November-March 

0.08 

0.08 

(per  kW  h) 

May-October  (per  kW  h) 

0.09 

0.09 

Daily  service  charge 

0.25 

0.25 

Gas  price  rate  (per  term) 

0.74 

0.74 

Monthly  service  charge 

9.70 

9.70 

The  use  of  nanofluids  in  solar  systems  is  desired  to  be  evaluated 
in  the  economic  point  of  view.  Otanicar  and  Golden  [26]  compared 
the  economic  features  of  nanofluid-based  solar  collectors  with  the 
conventional  types.  The  study  was  based  on  life  cycle  assessment 
which  is  a  capable  methodology  to  evaluate  the  economic  impacts 
of  products.  As  it  is  shown  in  Table  1,  the  capital  and  maintenance 
costs  are  $120  and  $20  higher  for  nanofluid-based  collectors 
respectively  for  the  year  life  time.  The  payback  period  is  less  for 
the  conventional  collectors,  but  with  assuming  15  years  life  time 
according  to  the  better  performance  of  nanofluid-based  collectors, 
the  life  cycle  savings  would  be  nearly  the  same.  Taylor  et  al.  [25] 
evaluated  the  utilization  of  nanofluid  receivers  for  a  solar  thermal 
power  plant  theoretically.  They  appraised  the  amount  of  nanopar¬ 
ticles  that  would  be  needed  for  a  solar  thermal  power  plant  (3  kg 
per  each  1  MWe)  and  with  assuming  the  price  of  nanoparticles 
~$1000/kg,  the  total  capital  investment  would  be  $5/W  where  the 
cost  rise  is  less  than  0.1%  of  the  total  capital  investments.  They 
proposed  two  conceptual  designs  for  nanofluid  receivers  that 
substitute  the  conventional  models  in  the  solar  plant  which  are 
cheaper  than  the  conventional  ceramic  receivers.  They  used 
graphite/therminol®  VP-1  nanofluid  in  the  model  with  a  volume 
fraction  less  than  0.001%.  The  conservative  calculations  exhibited 
that  for  a  100  MWe  power  tower  solar  plant  operating  in  Tucson- 
Arizona,  more  than  about  $3.5  million  was  obtained  in  the  yearly 
revenue  and  the  payback  time  of  the  plant  reduces  about  two 
years. 

The  objective  of  the  current  research  is  to  model  a  trough 
collector  absorber  tube  and  study  the  effects  of  nanoparticle 
concentration  on  the  mixed  convection  heat  transfer  rate  of  the 
nano-fluid  in  a  fully  developed  turbulent  flow.  Most  of  the  models, 
in  the  literature,  for  simulating  the  flow  in  a  trough  collector 
absorber  tube  assume  that  the  solar  flux  is  uniform  and  many 
correlations  in  the  models  are  based  on  a  uniform  temperature. 
Because  of  the  non-uniform  nature  of  the  solar  flux  on  the  outer 
surface  of  the  inner  absorber  tube,  the  temperature  distribution 
of  the  heat  transfer  fluid  is,  in  fact,  asymmetric.  Relatively  few 
investigations  have  been  reported  on  simulating  non-uniform  heat 
flux  in  absorber  tubes  [27,28].  Al203/synthetic  oil  was  used  as  the 
heat  transfer  fluid.  The  fluid  properties  are  strongly  dependent  on 


the  absorber  operational  temperature.  The  heat  transfer  coeffi¬ 
cients  for  various  volumetric  concentrations  of  nanoparticles  in 
the  heat  transfer  fluid  are  presented  along  the  axial  and  circum¬ 
ferential  directions  for  a  number  of  absorber  tube  operational 
temperatures. 


2.  Fluid  properties 

Thermo-physical  properties  of  the  base  fluid  (synthetic  oil)  and 
the  nanoparticles  (A1203)  including  density,  viscosity,  thermal 
conductivity,  and  specific  heat  are  strongly  dependent  on  the 
operational  temperature.  The  variation  of  density  of  the  base  fluid 
with  temperature  was  extracted  from  Ref.  [29]  and  plotted  in  Fig.  1(a). 
The  distributions  of  viscosity,  thermal  conductivity,  and  specific  heat 
of  the  base  fluid  as  a  function  of  temperature  were  also  obtained  from 
[27]  and  used  in  the  current  study  (Fig.  1(b)  and  (c)).  Also,  the 
properties  of  A1203  were  taken  from  Ref.  [30].  Addition  of  nanopar¬ 
ticles  remarkably  alters  the  thermo-physical  properties  of  the  heat 
transfer  fluid.  These  properties  were  calculated  for  various  concentra¬ 
tions  of  nanoparticles  in  the  fluid  and  plotted  in  the  same  figure  for 
compassion.  Table  Al  of  Appendix  A  provides  the  numerical  values  of 
these  properties  for  the  operational  temperatures  of  300  K,  400  K, 
and  500  K.  The  thermo-physical  properties  of  the  base  fluid,  nano¬ 
particles  and  nanofluid  as  a  function  of  temperature  are  provided  in 
Appendix  A. 


3.  Numerical  modeling 

The  flow  of  heat  transfer  fluid  through  an  LS-2  parabolic  trough 
collector  module,  tested  on  the  AZTRAK  rotating  test  platform  at 
SNL  by  Dudley  et  al.  [31],  was  numerically  simulated.  A  schematic 
diagram  of  the  collector  and  absorber  tube  is  shown  in  Fig.  2  and 
the  geometric  characteristics  of  the  collector  are  shown  in  Table  2. 
The  absorber  tube  is  7.8  m  (L)  long  and  0.07  m  ( Da )  in  diameter. 
Al203/synthetic  oil  was  used  as  the  working  fluid  in  the  collector 
field  and  the  nanofluid  flow  was  modeled  as  single  phase.  Since 
the  Re  is  higher  than  2500  and  Richardson  number  (Ri=Gr/Re2)  is 
between  0.1  and  10,  the  flow  is  assumed  to  be  fully  developed 
turbulent  mixed  convection.  Gravitational  acceleration  was  also 
taken  into  account.  The  operational  temperatures  are  300  K,  400  K, 
and  500  K  and  the  pressure  is  assumed  20  psig.  The  fluid  motion 
governing  equations  were  solved  using  the  FLUENT  software.  The 
governing  equations  are  provided  in  Appendix  A. 

A  k-e  formulation  was  used  to  model  the  turbulence  in 
the  flow.  In  the  k-e  model,  the  turbulence  characteristics  are 
described  based  on  the  turbulent  kinetic  energy  (TKE)  and  viscous 
dissipation  of  turbulent  kinetic  energy.  The  FLUENT  software  was 
used  to  numerically  simulate  the  flow  of  the  nano-fluid  in  the 
absorber  tube.  This  software  is  based  on  finite  volume  method. 
To  couple  the  velocity  field  and  pressure,  the  SIMPLE  model  was 
employed.  In  this  method,  the  pressure  correction  equation  is 
solved  in  a  number  of  iterations,  and  then  the  velocity  is  corrected 
until  the  continuity  equation  is  satisfied  in  the  computational 
domain.  For  the  convective  and  diffusive  terms,  a  second  order 
discretization  method  was  used.  The  flow  was  assumed  three- 
dimensional,  steady  state,  and  single-phase.  To  insure  that  the 
results  were  time  independent,  the  average  heat  transfer  coeffi¬ 
cient  was  monitored  with  the  size  of  the  grid. 

The  solar  energy  flux  distribution  on  the  outer  wall  of  the  inner 
absorber  tube  was  obtained  using  the  Monte  Carlo  Ray  Trace 
(MCRT)  technique  developed  by  He  et  al.  [27],  The  heat  flux 
distribution  of  He  et  al.  (Fig.  3)  was  used  in  the  current  study.  In 
the  numerical  simulations,  the  boundary  condition  on  the  absor¬ 
ber  tube  wall  was  set  to  this  non-uniform  heat  flux  distribution 
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T  (K)  T(K) 

Fig.  1.  (a)  Density,  (b)  viscosity,  (c)  thermal  conductivity,  and  (d)  specific  heat  of  the  base  fluid  as  well  as  the  nanofluid  for  various  particle  concentrations  at  different 
operational  temperatures. 


Fig.  2.  Schematic  diagram  of  the  parabolic  trough  collector  and  absorber  tube.  Red 
lines  represent  rays  of  light  from  the  Sun.  (For  interpretation  of  the  references  to 
color  in  this  figure  caption,  the  reader  is  referred  to  the  web  version  of  this  article.) 


e 


Fig.  3.  Heat  flux  distribution  of  the  absorber  tube  as  a  function  of  angle.  Graph  was 
reproduced  from  [27]. 


Table  2 

Specifications  of  the  LS-2  parabolic  trough  collector  [31]. 


System  Type 

Parameter 

Value 

Parabolic  trough  reflector 

Z,*x, 

7.8  m*5  m 

f 

1.84  m 

Receiver 

h 

0.0575 

ra 

0.035 

using  a  user-defined  function  in  the  FLUENT  software.  In  addition, 
the  boundary  conditions  at  the  inlet  (z= 0)  were  as 

w  =  w0,  u  =  v  =  0,  T  =  Tj 

where  u,  v ,  and  w  are  the  velocity  components  in  x,  y,  and  z 
directions,  respectively,  and  T  is  the  temperature  of  the  fluid.  The 
flow  at  the  outlet  of  the  tube  was  assumed  fully  developed  and  an 
outflow  condition  was  used  in  the  simulations,  the  fully  developed 


- Z/D=2 

- Z/D=8 

- Z/D=30 

.  Z/D=60 

a  Z/D=100 
-  Z/D=l  10 


Fig.  4.  Axial  velocity  profiles  at  different  axial  positions  at  Re=26  000. 


entrance  length  is  demonstrated  in  Appendix  A.  Also,  Fig.  4  shows 
the  axial  velocity  profiles  at  six  different  cross  sections  for  a  given 
Re  (260  00).  We  can  observe  that  the  fully  developed  region  occurs 
in  Z/D =30. 
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Furthermore,  the  following  boundary  conditions  were  used  at 
the  fluid-solid  interface: 

n  »  dT 

u  =  v  =  w  =  0,  q"=-Keff— 

The  MCRT  method  is  dependent  on  the  geometric  (r,  ZL )  as  well 
as  the  grid  ( Nc ,  Nz )  parameters.  The  MCRT  method  requires  a  grid 
density  for  the  absorber  tube  of  Nc  >  68  and  Nz  >  320.  A  structured 
non-uniform  mesh  was  created  to  discretize  the  computational 
domain.  It  is  finer  near  the  wall  where  the  velocity  and  tempera¬ 
ture  gradients  are  larger.  The  standard  method  to  test  for  grid 
independence  is  to  increase  the  resolution  and  repeat  the  simula¬ 
tion.  If  the  results  do  not  change  appreciably,  the  original  grid  is 
probably  adequate  and  the  results  are  not  dependent  on  the  girds 
size.  If,  on  the  other  hand,  there  are  significant  differences 
between  the  two  solutions,  the  original  grid  is  likely  of  inadequate 
resolution.  In  such  a  case,  an  even  finer  grid  should  be  tried  until 
the  grid  is  adequately  resolved.  A  number  of  grid  sizes  were 
studied  to  insure  that  the  accuracy  of  the  results  is  independent  of 
the  grid  size.  The  mean  heat  transfer  coefficient  was  monitored  for 
two  cases:  (1)  T=500K  and  cp= 0  and  (2)  T=500K  and  q>=0.05. 
The  results  of  the  grid  refinement  study  are  presented  in  Table  3. 
Increasing  the  grid  numbers  does  not  significantly  change  the  mean 
heat  transfer  coefficient.  The  final  mesh  for  the  current  study 
consisted  of  68  (Nc),  320  (Nz),  and  20  (Nr)  nodes  in  the  circumfer¬ 
ential,  axial  and  radial  directions,  respectively. 


Table  3 

Average  heat  transfer  coefficients  at  500  K  for  various  grid  sizes. 


Grid  (l\l?N?l\y 

Mean  heat  transfer 
coefficient  at  <p=0 

Mean  heat  transfer 
coefficient  at  </>=0.05 

68-20-320 

137.99 

146.53 

68-30-400 

137.75 

147.01 

80-20-480 

138.01 

146.55 

4.  Validation 

Two  available  experimental  cases  were  used  to  validate  the 
current  numerical  study. 

4.2.  LS-2  solar  collector 

Three  typical  states  for  a  bare  tube  from  Dudley  et  al.'s  report 
[31]  were  used  for  the  purpose  of  validation.  The  results  are 
presented  in  Table  4.  Dudley  et  al.  derived  the  performance 
correlation  between  the  collector  efficiency  and  the  working  fluid 
temperature  as  following  [31]: 

r\  =  I<*( 73.3  -  0.007276*AT)+  ^4.96*^ j  j  -  0.0691* 

In  the  equation,  I<  is  the  incident  angle  modifier  and  may  be 
assumed  to  be  one  (2<  =  1). 

The  heat  flux  distribution  on  the  outer  surface  of  absorber  tube 
is  shown  in  Fig.  5(a).  The  temperature  distribution  on  the  outer 
surface  of  absorber  tube  is  also  shown  in  Fig.  5(b).  Fig.  5(b)  shows 
the  wall  temperature  distributions  of  the  absorber  tube  based  on 
the  test  condition  state  1.  Fig.  6  shows  the  temperature  distribu¬ 
tion  in  the  axial  direction  with  6=0°,  9CT,  180°,  270°  and  310°.  The 
wall  temperature  increases  with  the  increase  of  z  along  axial 
direction,  and  the  temperature  at  310°  is  much  larger  than  the 
other  angular  direction.  The  reason  is  that,  0=310°  is  on  the  side 
with  high  concentrated  solar  radiation,  while  0=90°  is  on  the  side 
without  concentrated  solar  radiation.  According  to  the  simulated 
outlet  temperature  and  collector  efficiency  based  on  the  test 
results  of  Dudley  (Table  4),  the  absolute  error  is  less  than  3.8  °C 
which  has  good  agreement  with  the  obtained  simulation  results. 
This  proves  that  the  models  and  methods  used  in  the  present 
study  are  feasible  and  the  numerical  results  are  reliable 


Table  4 

The  experimental  data  from  Dudley  and  comparisons  with  our  simulation  result. 


Test 

4  (W/m2) 

m  (kg/s) 

Reynolds 

Ta  (”C) 

Tin  (°C) 

To.test  (°C) 

To.model  (°C) 

£(%) 

4(test)  (^0 

4 (model)  (^) 

State  1 

933.7 

0.66 

4501 

21.2 

102.2 

124 

126.1 

1.69 

72.51 

72.27 

State  2 

968.2 

0.65 

7653 

22.4 

151 

173.3 

176.9 

2.1 

70.9 

71.5 

State  3 

937.9 

0.61 

28560 

28.8 

297.8 

316.9 

319.32 

0.76 

68.98 

66.9 

Fig.  5.  (a)  Heat  flux  and  (b)  temperature  distributions  on  the  outer  surface  of  the  absorber  tube. 
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4.2.  Al203  nanofluid  flow  through  a  horizontal  tube 

Heat  transfer  coefficient  and  Nusselt  number  were  obtained  and 
compared  with  the  experimental  results  of  Esmaeilzadeh  et  al.  [32] 
in  a  horizontal  tube.  Fig.  7(a)  and  (b)  shows  the  both  the  results  of 
the  current  numerical  study  and  those  of  Esmaeilzadeh  et  al.  [32], 
The  axial  evolution  of  Nusselt  number  and  heat  transfer  coefficient 
are  in  good  agreement  with  the  corresponding  experimental  results. 


z 


Fig.  6.  Mean  temperature  distributions  along  the  axial  direction  (z)  on  the  outer 
surface  of  absorber  tube  at  various  angular  locations  at  Re =17  361. 


5.  Results  and  discussion 

The  study  of  heat  transfer  was  carried  out  on  a  trough  collector 
absorber  tube  with  7.8  m  length  and  0.07  m  diameter.  The  tube  is 
consumed  to  be  evacuated  inside  the  space  between  the  absorber 
metal  tube  and  the  glass  tube.  The  vacuum  level  in  the  annular 
space  between  the  inner  absorber  tube  and  the  outer  glass  cover 
tube  is  so  high  that  the  convective  heat  losses  could  be  essentially 
eliminated,  but  the  thermal  radiation  heat  transfer  between  the 
outer  wall  of  absorber  tube  and  the  inner  surface  of  glass  cover 
tube  is  taken  into  account.  Also,  the  heat  transfer  mechanism 
inside  the  tube  is  as  convection  and  the  inlet  fluid  velocity  is 
assumed  to  be  4  m/s.  The  convection  heat  transfer  of  nanofluid 
may  be  increased  compare  to  the  base  fluid.  The  improved  thermal 
conductivity  of  nanofluid  is  the  most  important  factor  for  improv¬ 
ing  heat  transfer  in  parabolic  trough  collectors.  Actually,  nanofluid 
reduces  the  thermal  resistance  at  interfaces  and  minimizes  the 
temperature  difference  between  the  absorber  and  the  heat  trans¬ 
fer  fluid.  The  effects  of  the  particle  concentration  on  the  flow  fields 
and  heat  transfer  coefficient  along  the  absorber  tube  are  pre¬ 
sented.  Convection  heat  transfer  coefficients  at  various  axial 
locations  for  a  given  Reynolds  number  and  five  volumetric  con¬ 
centrations  of  A1203  particles  are  shown  in  Fig.  8.  Because  of  the 
non-uniform  heat  flux  nature,  a  number  of  graphs  are  shown  for 
heat  transfer  coefficients  (Fig.  8).  As  shown,  the  heat  transfer 
coefficient  ranges  from  100  to  260  and  60  to  160  at  the  operational 
temperatures  of  300  K  and  500  K,  respectively,  indicating  the 


Fig.  7.  Comparison  of  (a)  the  heat  transfer  coefficients,  (b)  the  Nusselt  number  along  the  axial  direction  from  the  current  numerical  study  with  the  existing  experimental 
values. 


Fig.  8.  Convective  heat  transfer  coefficient  vs.  axial  distance  for  various  particle  concentrations  and  angular  locations  for  (a)  Re=3900  (T=300K)  and  (b)  Re=26  800 
(1=500  K). 
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Fig.  9.  Mean  convective  heat  transfer  coefficient  vs.  particle  concentration  at  the 
operational  temperatures  of  300,  400  and  500  K. 


T  (K) 

Fig.  10.  Percentage  increase  in  mean  convection  heat  transfer  coefficient  due  to  the 
addition  of  nanoparticles  to  the  base  fluid  vs.  operational  temperature  for  various 
particle  volume  concentrations. 


decrease  of  heat  transfer  coefficient  with  the  increase  in  opera¬ 
tional  temperature. 

The  heat  transfer  coefficient  sharply  drops  and  then  levels 
off  to  its  fully  developed  value.  As  the  nanoparticle  concentration 
increases,  the  heat  transfer  coefficient  rises.  The  coefficient  is  more 
affected  by  the  particle  concentration  at  the  operation  tempera¬ 
ture  of  300  I<  than  500  K.  However,  the  heat  transfer  coefficient  is 
reduced  as  the  operational  temperature  is  increased  from  300  K  to 
500  K.  The  average  heat  transfer  coefficients  are  lower  for  the 
higher  temperatures  (Fig.  9).  As  shown  in  this  figure,  the  average 
convection  heat  transfer  coefficient  increases  by  about  14%  by 
increasing  the  nanoparticle  volumetric  concentration  by  5%  at  the 
operational  temperature  of  300  K.  Also,  the  effect  of  nanofluid 
concentrations  on  average  heat  transfer  coefficient  at  the  opera¬ 
tional  temperature  of  400  K  increases  by  8.6%  with  the  addition  of 
5%  (in  volume)  of  A1203  nanoparticles.  This  effect  was  also  studied 
at  500  K  and  the  results  show  that  the  average  convective  heat 
transfer  coefficient  increases  6%  by  adding  5%  (in  volume)  of  A1203 
nanoparticles. 

From  Fig.  1,  we  can  observe  that  when  the  operational 
temperature  is  300  K  the  value  of  / 1  increases  to  8.38e-4  and 
the  value  of  k  increase  to  0.028  by  adding  5  vol%  of  A1203 
nanoparticles.  While  at  the  operational  temperature  of  400  K, 
the  value  of  /(  increases  to  2.7e-4  and  the  value  of  k  increase  to 
0.023  by  adding  5  vol%  of  A1203  nanoparticles.  Also  at  the  opera¬ 
tional  temperature  of  500  K  the  value  of  increases  to  9.6e  -  5 
and  the  value  of  k  increases  to  0.019  by  adding  5  vol%  of  A1203 
nanoparticles.  The  other  reason  for  the  lower  increase  in  heat 
transfer  coefficient  at  higher  temperatures  can  be  associated  to  the 


decrease  in  viscosity  and  thermal  conductivity  change  rates  at 
higher  temperatures  [23],  The  percentage  increase  in  heat  transfer 
at  various  temperatures  is  shown  in  Fig.  10.  According  to  this 
figure,  the  percentage  of  heat  transfer  increase  has  a  descending 
trend  with  temperature  rising. 


6.  Conclusions 

The  thermal  behavior  of  a  nanofluid  consisting  of  synthetic  oil 
and  A1203  throughout  a  parabolic  trough  collector  absorber  tube 
(L= 7.8  and  D= 0.07  m)  has  been  simulated  and  a  turbulent  mixed 
convection  heat  transfer  with  non-uniform  heat  flux  has  been 
investigated.  The  results  of  the  current  study  illustrate  that  the 
heat  transfer  coefficient  of  the  working  fluid  in  an  absorber  tube  is 
enhanced  at  the  presence  of  nanoparticles.  The  heat  transfer 
coefficient  is  increased  as  the  concentration  of  the  nanoparticles 
in  the  base  fluid  is  increased.  For  a  given  inlet  Reynolds  number, 
the  enhancement  in  heat  transfer  coefficient  caused  by  nanopar¬ 
ticles  is  decreased  as  the  operational  temperature  of  the  absorber 
tube  is  increased.  The  maximum  heat  transfer  coefficients  were 
obtained  at  the  left  and  right  sides  of  the  tube  (0=310  and 
0=23O).The  results  of  this  work  for  investigating  the  heat  transfer 
variations  of  the  assumed  collector  were  validated  with  the  other 
similar  studies. 

Utilizing  nanofluid  in  PTC  has  environmental  benefit  and 
causes  improving  heat  transfer,  consequently  reducing  the  needed 
heat  transfer  area  of  the  tubes  and  heat  exchangers.  Nevertheless, 
these  advantages  could  be  achieved  with  little  change  (in  terms  of 
materials,  system  design,  and  initial  capital  investment)  on  the 
entire  PTC  system. 


Appendix  A 


Fluid  properties: 

pbf  =  -0.9985*1+1236  kg/m3 

(1) 

ps  =  3850  kg/m3 

(2) 

Peff  =  ^-tP)Pbf  +  <PPs 

(3) 

The  correlation  of  effective  density  is  taken  of  [33]: 

Py  =  6.67*10-7T4  — 1.56*1  CP3T3  +  1.38*T2  —  5. 541  *102r+ 8. 848*1 04  /(Pa  s 

(4) 


=  (5) 
The  correlation  of  effective  viscosity  is  taken  of  Einstein's  [34]: 
Kbf=  -5.753496*10-10r2-1.875266*10-4T  +  1.900210*10-'  W/mK 

(6) 

Ks  =  5 . 5  +  34. 5*e(  ~  0  0033*Cr  -  273»  w/m  K  (7) 


Keff  =  Kf* 


Ks  +  2Kf-2(\+fi)3(p{Kf-Ks) 
Ks  +  2Kf  +  (\  -Ks) 


(8) 


where  fi  is  the  ratio  of  nanolayer  thickness  to  the  nanoparticles 
diameter  which  is  considered  to  be  0.1  in  Eq.  (8)  of  this  study  [35], 


Cpbf  =  0.0017081+1.107798  kj/kg  K 


(9) 


Cps  =  1.046+1.74*10~4r-2.79*104r2  kj/kg I<  (10) 


(Cp)eff  =  a-<P)Cpbf  +  tpCPs 


(ID 


The  correlation  of  effective  specific  heat  is  taken  of  [33], 
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Table  A1 

Nanofluid  (A1203)  thermo-physical  properties  at  300,  400  and  500  K. 


Property 

7=300  K 

7=400  K 

7=500  K 

0 

1% 

3% 

5% 

0 

1% 

3% 

5% 

0 

1% 

3% 

5% 

P  (kg/m3) 

936.5 

965.6 

1023.9 

1082.1 

836.6 

866.7 

927.0 

987.3 

736.7 

766.9 

830.1 

892.4 

Cp  0/kg  I<) 

1620 

1612 

1595 

1578 

1791 

1782 

1765 

1748 

1962 

1952 

1934 

1915 

K  (W/m  K) 

0.134 

0.139 

0.150 

0.161 

0.115 

0.119 

0.129 

0.138 

0.096 

0.100 

0.107 

0.115 

pi  (mPa  s) 

6.68 

6.84 

7.18 

7.518 

2.16 

2.22 

2.33 

2.43 

0.77 

0.789 

0.828 

0.866 

In  all  of  the  correlations  /,  s  and  eff  correspond  to  the  based 
fluid,  nanoparticle  and  nanofluid,  respectively  and  (p  is  the  volu¬ 
metric  concentration  of  A1203  in  the  nano-fluid.  The  thermo¬ 
physical  properties  of  the  nanofluid  were  calculated  via  (3),  (5), 
(8),  (11).  The  governing  equations  for  the  fluid  motion  in  a 
cylindrical  coordinate  system  are  described  as  follows: 

Continuity  equation: 

7^/V/u>+7^(/V/rv>+^%/W)  =  0  (12) 

Momentum  equation: 

8  component: 

\  j0(PeffUU)  +  l  'l(PeffrVU)  +  ^(Peffwu)  +  \peffUV) 

1  fdp\  1  d  ("  du\  '  d  (heffd(ru)\ 

~  r\ddj+r 2  d 8  ye)id8 J  dr\  r  dr  J 

■+ — ^r2f  dQ+PeffSPeff^  W  -  T)  sin  8  (13) 

r  component: 

\td{Pnfuv)+\i{Pnfrvv)+tz(PnfWV)  ~  l(P"fuI) 

l/dp\  1  d  /  <5v\  d  (P„fd(rv)\ 

~  r\dr J  +  r2  d0\Mnf  dd J  +  dry  r  dr  J 

-  ^  -  PnfSfi„f(J w~T)  cos  0 

2  component: 

7^e//uw)+l^(pe//rvw)+^(pe//w2) 

dp  Id/  dw\  Id/  dw\ 

=  ~^+r2W\Peffl)0)  +  rdr\rPeff^dr ) 

Energy  equation: 

ld0(peffuT)  +  ldr(pe»rvT) + l^P'flwT) 

=  1A(  Keff  dT\  d  (  Keff  dT\ 
r2  d8  \(Cp)eff  d0J  rdr  ^  (Cp)eff  dr  J 

where  r,  8,  and  z  are  the  coordinates,  u,  v,  and  w  are  the  velocity 
components  and  p,  T  and  g  are  the  pressure,  temperature  and 
gravitational  acceleration,  p,  k,  p,  and  Cp  are  the  effective  density, 
thermal  conductivity,  dynamic  viscosity  and  specific  heat  capacity 
of  the  nano-fluid. 

The  fully  developed  entrance  length: 

1.22  <Lh  =  (4.4*Re’ /6  )*D  <  1 . 

6Lt=  10*D  =  0.7  m 

See  Table  Al. 
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